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The supercurrent of a Josephson junction is reduced by phase diffusion. For ultrasmall capaci- 
tance junctions the current may be further decreased by Coulomb blockade effects. We calculate the 
Cooper pair current by means of time-dependent perturbation theory to all orders in the Josephson 
coupling energy and obtain the current-voltage characteristic in closed form in a range of parameters 
of experimental interest. The results comprehend phase diffusion of the coherent Josephson current 
in the classical regime as well as the supercurrent peak due to incoherent Cooper pair tunneling in 
the strong Coulomb blockade regime. 
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New lithography and low-temperature techniques have 
allowed the fabrication and measurement of small 
Josephson junctions affected by the capacitive charging 
energy of single Cooper pairs [[[j-p). Much of the work 
so far has concentrated on the region of strong Coulomb 
blockade where the tunneling of Cooper pairs described 
by the Josephson energy can be treated perturbatively. 
While Coulomb blockade of Cooper pair tunneling is 
fairly well understood, the relation between the effects 
observed at low temperatures and the familiar "classical" 
dynamics of Josephson junctions B remains to be exem- 
plified. In this article we demonstrate that the current 
peak caused by incoherent Cooper pair tunneling in the 
regime of strong Coulomb blockade gradually evolves into 
the classical supercurrent when parameters are changed 
accordingly. 

We consider a Josephson junction with capacitance C 
and critical current I c = (2e/fr)Ej, where Ej is the 
Josephson energy. The capacitance gives rise to a charg- 
ing energy E c — 2e 2 /C for Cooper pairs. The junction 
is coupled to an ideal voltage source through a resistor 
of resistance R as shown in Fig. [p. This system can be 
modeled by the Hamiltonian 



H — Hj + H env 
where the first term 



H, 



-Ej cos(». 



(1) 



(2) 



is related to the tunneling of Cooper pairs through the 
Josephson junction. In analogy to the familiar relation 
between the voltage across the Josephson junction and 
the corresponding phase difference 



we may introduce a phase 



ip R = —Vt - <p 



(4) 



related to the voltage Vr = V — Vj across the resistor 
where V is the external voltage. In the following, we are 
interested in a voltage regime far below the gap voltage 
so that tunneling of quasiparticles may be neglected. 

The second part of the Hamiltonian ([j]) describing the 
coupling of the junction to its electrodynamic environ- 
ment reads H 



H 



Q 2 



2C ^ 



2C n 



(5) 



The first term represents the charging energy where the 
charge Q obeys the commutation relation [ip, Q] = 2ie. 
By virtue of Eq. (Eh , the second part describes a bilinear 




(3) 



FIG. 1. A Josephson junction characterized by the Joseph- 
son energy Ej and capacitance C is coupled to a voltage 
source V via a resistor R. The voltage drops across the junc- 
tion and the resistor are Vj and Vr, respectively. 
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coupling between the phase difference across the Joseph- 
son junction and phases of a set of LC-circuits which 
model the electromagnetic environment of the junction. 
The phases ip n and charges q n are conjugate variables 
according to [ip n ,q n ] — 2ie. It can be shown that by 
suitably choosing the parameters C n and L n , the second 
part of H cnv effectively describes the resistor R. 

The thermal equilibrium Cooper pair current at tem- 
perature T — l/fcs/3 may be written as 



I c (W + (oo, t ) sin[^(t )Moo, to)) 
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(6) 



where the time dependence of ip(t) arises from the Hamil- 
tonian (0). 



U(t,t ) = T exp (^J^ dt'Ejcos[ip(t')] 



(7) 



with the time ordering operator T is the time evolution 
operator in the interaction representation. Taking the 
limit t — > oo, the precise value of the initial time to be- 
comes irrelevant and will be set to zero in the following. 
Expanding the current (||) in powers of the Josephson 
energy Ej one obtains 



* -'IS 

M=l 



2h 



-E., 



2M-1 




2iW-l 



n ^ ^ 



fe=i 



x / dti 
'o 



dt 



2M-1 



,2e 



'0 

2M-1 



exp i— V ^2 Cfeife 



x exp 



k=0 
2M-1 fc-1 

EE 

k=l 1=0 



CkClJ[Vk(tk -ti)] 



(8) 



For given order M, the Cfc, k = 0, . . . , 2M — 1, which arise 
from a decomposition of the trigonometric functions in 
Eqs. (H) and (Q) into exponentials, may take the values 

±1 and are subject to the constraint Y^=o 1 ^ k = ®- ^ ne 
factors r\ k also take the values ±1 and account for the 
fact that Eq. (|^) contains two time-ordered operators. 
Finally, the phase correlation function || 



J(t) 



du> KeZf (lu) e 



R Q 1 - 



(9) 



contains the relevant information on the total environ- 
mental impedance. 

In our case, the impedance as seen from the tunnel 
junction contains a capacitance C in parallel with a re- 
sistance R leading to 



ReZ t (u) 
Rq 



1 + (lu/lorY 



(10) 



The dimensionless impedance at zero frequency is p — 
R/Rq with the resistance quantum Rq = h/4e 2 , while 
the frequency scale is set by lor = 1/RC. 

For the impedance ( |l0| ) , evaluation of the integral 
gives 



J(t) = -2p 



E 



! e -"R\t\ 



-Vn\t\ 



cot 



+.2(1 



^ n[l - (j/„/wr) 2 ] 
-^l*l)sgn(t) 



with 



5 = 7 + 



2 

7T (0 



(11) 



f!2) 



Here, 7 = 0.5772 ... is Euler's constant and the v n — 
2Trn/Ti(3 are the Matsubara frequencies. Further, ip(x) 
is the logarithmic derivative of the gamma function and 
sgn(x) denotes the signum function. 

Now, typical lead resistances are of the order of 100f2 
H so that p I. Then, for small Josephson junctions 
the large lead conductance overdamps the junction and 
ujr is large compared to the Josephson frequency 



2e 



-RIc- 



(13) 



Hence, we may disregard terms proportional to 
exp(— ujRt), and Eq. (11) reduces to 



J(t) = -2p 



hp 



\t\ + s-j2 



-v n \t\ 



-fi-sgn(i) 



(14) 



Furthermore, since p <C 1, the ratio v\/loj — l/pf3Ej 
may be large even for low temperatures, and the result 
( |l4| ) can be simplified further to read 



J(t) = -2p 



fyt\ + S + q sg n(t) 



(15) 



Inserting this correlation function into (g), the time in- 
tegrals may be performed and one obtains the expansion 



a 2M-1 

Is = / c Re [ J2 E II a ( x *) ] > 

M=l{x k } k=l 



(16) 



where the set of positive integers {x k } describes a se- 
quence with \xk+i — x k \ = 1 starting and ending at 
xi = X2M—1 = !• Each value of Xk is associated with 
a factor 



a(x k ) 



sin(irpx k ) exp[2pSx k (x k+ i - x k )] 



2npx k 



v + ix k //3Ej 



(17) 
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which contains the dependence of the Cooper pair current 
on the applied voltage V via the dimensionless variable 
v = V/RI C . An expansion of the form ( |l6| ) appears in 
the context of the one-dimensional Coulomb gas and 
was also found in a previous treatment of the dynam- 
ics of small current-biased Josephson junctions || for a 
different kind of approximation. Following H], we may 
rewrite ( |l6| ) in terms of a continued fraction 



Is = I c Rc 



sin(7rp) cxp(— 2pS) 



2np v + i/pEj 1 



61 



(18) 



with coefficients 



' (3Ej\ sin(7rpn) sin(7rp(rc. + 1)) exp(— 2pS) 
1 2^p~) n(n + l)(n - w(3Ej)(n + 1 - iv(3E,j) ' 



(19) 



It is instructive to first investigate the limit p 
where the coefficients simplify to read 



b n = 



(n - iv(3Ej){n + 1 - iv/3Ej) 



(20) 



and the continued fraction (jig ) may be evaluated by 
means of a matrix recursion S. The recursion relations 
can be solved in terms of modified Bessel functions of 
complex order. For finite p <C 1 the corrections are not 
necessarily small since in the current-voltage characteris- 
tic (|l^) p is multiplied by a factor S defined in Eq. Jl^ ) 
which may contain large terms of order \og((3E c / p). On 
the other hand, a numerical evaluation of the continued 
fraction converges rapidly and only coefficients b n with 
moderately large n are relevant. Thus, for p <C 1, the 
sine functions in Eq. ( |l9|) may be linearized and it is 
readily seen that the dominant corrections are described 
in terms of an effective Josephson energy 



E*j = EjeM-pS)- 



(21) 



The continued fraction may then be evaluated by the 
method described above leading to the current-voltage 
characteristic 



Il-if3eV/7T P {f3Ej) 
I-if3eV/Trp((3Ej) 



(22) 



This constitutes the central result of this paper. The 
dependence of the effective Josephson energy on temper- 
ature is depicted in Fig. 2 for various values of p. 

To investigate the classical limit of the current-voltage 
characteristic (E2|) one must know that the limit h — > 
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f3E c 



FIG. 2. The effective Josephson energy Ej is shown as a 
function of the inverse temperature for p = 0.01, 0.05, and 0.1 
from the upper to the lower curve. 

has to be performed such that the flux quantum h/2e re- 
mains constant. This means that p is of order Ti and E c of 
order Ti 2 . Then, in the overdamped limit u> r 3> wj the ef- 
fective Josephson energy coincides with the bare Joseph- 
son energy Ej. This can also be seen in Fig. || where for 
small p the range in which Ej practically equals Ej be- 
comes very large. In this case, the current-voltage char- 
acteristic (22) can be shown to reduce to the standard 
Ivanchenko-ZiPberman result |l(],[ll]] for classical over- 
damped Josephson junctions. 

While in the classical limit (3E C / p vanishes, in the limit 
of strong Coulomb blockade we have (3E C / p ^> 1. The 
psi function in Eq. ( |l2] ) may then be approximated by a 
logarithm to obtain 



S = 1 + \n{{3E c /2ir 2 p). 



(23) 



Current-voltage characteristics for overdamped ultra- 
small Josephson junctions have so far only been calcu- 
lated to lowest order in the Josephson energy. The su- 
percurrent arising from incoherent Cooper pair tunneling 
then reads 



/ 



\T{p-t(3eV/^ 
T{2p) 



sinh(/3eV) 



where 



J h E C P \ 2tt 2 



l-2p 



exp(-2p7). 



(24) 



(25) 



In Fig. ||, we compare our result (^2|) (full line) with 
the standard Ivanchenko-Zil'berman result (dotted line) 
and the Coulomb blockade result (|24|) (dashed line) for 
(3Ej = 1, p = 0.04, and two different values for (3E C . 
For large /3E C , the current- voltage characteristic ( p^ ) 
is in very good agreement with the Coulomb blockade 
result (24). With decreasing f3E c , the results (|22|) and 
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^1 




V/RI C 

FIG. 3. The current-voltage characteristic of a Josephson 
junction with Josephson energy (3Ej = 2 is shown for charg- 
ing energies {3E C = 1 and 10 5 and external resistance p = 0.04. 
The full line corresponds to our result ( |2^ ) while the dotted 
line gives the standard Ivanchenko-ZiPberman result and the 
dashed line depicts the prediction (^ij) for strong Coulomb 
blockade. The two current-voltage characteristics are verti- 
cally shifted with respect to each other by I/I c = 0.3 for sake 
of clarity. 



son coupling energy Ej. Noting that ultrasmall junc- 
tions in a standard electromagnetic environment are over- 
damped, i.e., loj/lo r = (2e/fi)R 2 I c C <C 1, and with the 
assumption p/3Ej = (eRI c /nkBT) <C 1, the general re- 
sult was written as continued fraction, which for p « 1 
could be summed in closed form. The inferred main re- 
sult (|2^) was shown to describe the changeover from the 
classical Josephson effect in the presence of phase diffu- 
sion to the quantum regime of strong Coulomb block- 
ade where Cooper pairs tunnel incoherently. The the- 
ory covers the experimentally relevant range for small 
capacitance Josephson junctions. Although detailed ex- 
perimental studies of the region between the above- 
mentioned limits are absent, recent work jl4| indicates 
quantum effects in qualitative accord with the predic- 
tions made. 

We would like to thank the authors of Ref. jl4| for 
stimulating discussions on experimental aspects of the 
problem and acknowledge helpful conversations with 
G. Goppert. Financial support was provided by the 
Deutsche Forschungsgemeinschaft . 



( p4| ) evolve differently and a crossover to the Ivanchenko- 
ZiPberman result can be found. 

To make the connection between our result ( |22| ) and 
the Coulomb blockade result ( p4] ) more explicit, we con- 
sider the current- voltage characteristic ( |22| ) to order E 2 
using the approximation (p3f). We then obtain for the 
Cooper pair current 



Is = f 



2n 2 p(3eV 
((3eV) 2 + ttV ' 



(26) 



While the approximation of the phase correlation func- 
tion J(t) leading from Eq. (|l4| ) to Eq. ( |l5|) is only justified 
for temperatures with p(3Ej <C I, for Ej <C E c the the- 
ory should extend to rather low temperatures. In fact, 
the zero bias differential resistance as well as the voltage 
for which the current takes its maximum are obtained 
precisely for small p including terms of order p 2 . 

The result ( p6|) may be improved by going beyond 
the approximation ( |l5|) for the phase correlation func- 
tion J(t). Keeping the sum over the Matsubara frequen- 
cies in Eq. (|l4|), a closed form for J(i) can be found in 
the strong Coulomb blockade regime (3E C / p ^ 1. The 
current-voltage characteristic (||) may then be evaluated 
to order E 2 and reproduces exactly the expression (|24|). 

In summary, we have derived a formally exact expan- 
sion of the current-voltage characteristic of a voltage- 
biased Josephson junction as a power series in the Joseph- 
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